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EXECUTIVE  SUMMARY 


An  account  is  provided  of  the  technical  work  carried  out  during  the 
period  between  delivery  of  the  preliminary  Solar  Wind  Transport 
Magnetosheath  Forecast  Model  (Phase  I  SWT  model)  in  November  1991  and 
the  early  termination  of  the  contract  in  July  1993. 

The  technical  work  carried  out  during  this  period  focused  on 
validation  of  the  SWT  model  and  has  involved  two  primary  technical 
tasks:  (1)  creation  of  the  observational  satellite  data  base  against 
which  the  model  was  to  be  validated,  and  (2)  initial  comparative  testing 
of  the  Phase  I  SWT  model  forecasts  against  these  observations.  Details 
of  the  results  obtained  from  these  two  technical  tasks  are  discussed. 

At  the  present  stage  of  development,  the  SWT  forecast  model  remains 
unqualified  for  operational  use.  An  outline  is  provided  of  the 
technical  tasks  required  to  complete  the  qualification  program. 
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1.  OVERVIEW  OF  SOLAR  WIND  TRANSPORT  MAGNETOSHEATH  FORECAST  MODEL 
1.1  Overall  Objectives 


The  overal'  objectives  of  the  current  contract  effort  Is  to  develop 
a  Solar  Wind  transport  Magnetosheath  Forecast  Model  (SWT  Model)  for 
operational  use  by  the  Air  Weather  Service  (AWS)  and  the  Space  Forecast 
Center  (SFC).  This  tnodel  will  be  the  first  component  In  a  sequential 
predictive  chain  that  will  comprise  an  Integrated  space  plasma  weather 
forecast  system  proceeding  from  Li  down  to  low  Earth  orbit. 

The  SWT  model  will  proceed  from  a  knowledge  of: 

•  real-time  solar  wind  properties  at  an  upstream  spacecraft 
monitor  located  approximately  at  the  Lj  llbratlon  point 


•  terrestrial  geomagnetic  field  and  certain  magnetospherlc 
Indices 


to  forecast  the  time-dependent  plasma  and  magnetic  field  properties 


•  from  the  upstream  monitor  to  and  through  the  terrestrial 
bow  shock 

•  throughout  the  3*0  volume  of  the  magnetosheath  region  down 
to  and  across  the  magnetopause  boundary 

•  Into  the  outer  region  of  the  magnetosphere 


The  operational  domain  of  the  Solar  Wind  Transport  Magnetosheath 
Forecast  model  Is  as  shown  in  the  sketch  below. 


.  J 

Operational  Domain  of 
Solar  Wind  Transport 
Magnetosheath  Forecast 
Model 
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Illustration  of  Operational  Domain  of 
Magnetosheath  Forecast  Model 


Solar  Wind  Transport 
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The  specific  output  from  the  SWT  model  will  be: 


•  location  of  bow  shock  and  magnetopause  surfaces 

•  plasma  and  magnetic  field  properties  throughout  the 
magnetosheath  region,  down  to  the  magnetopause,  and  into 
the  outer  magnetospherlc  region 

•  time  histories  of  plasma  and  magnetic  field  properties 
along  specific  spacecraft  trajectories 

•  input  data  for  Hagnetospheri c  Specification  Model  and 
Magnetopause  Boundary  Layer  Model 

•  selected  output  displays  to  facilitate  operational  use 


1.2  Operational  Elements  of  SWT  Model 

The  primary  challenge  in  the  development  of  the  Solar  Wind 
Transport  Magnetosheath  Forecast  Model  is  to  develop  an  operational 
forecasting  capability  to  be  approximately  100  times  faster  than  the 
actual  events  within  the  operational  domain  of  the  SWT  model  occur.  This 
will  require  SWT  model  operational  times  of  less  than  a  minute,  and  will 
provide  overall  warning  times  from  approximately  10  to  50  minutes 
depending  upon  oncoming  solar  wind  conditions. 

The  objective  of  the  first  18  months  of  the  planned  five-year 
contractual  effort  was  to  initially  establish  all  of  the  essential 
operational  elements  of  Solar  Wind  Transport  Magnetosheath  Forecast 
Model  and  to  embody  them  Into  a  preliminary  SWT  model.  These  elements 
are : 


1.  Input  Element 

•  measurement  of  solar  wind  properties  by  upstream  monitor 
located  at  the  Lj  libration  point 

•  determination  of  input  parameters  from  monitor  data  for 
solar  wind  transport  forecast  model 


2.  Solar  Wind  Transport  Forecast  Element 

•  determination  of  plasma  and  magnetic  field  properties 
between  upstream  monitor  and  bow  shock 

•  determination  of  forecast  time  for  Magnetosheath  and  Bow 
Shock  Forecast  Element 


•  determination  of  Input  parameters  for  Magnetosheath  and 
Bow  Shock  Forecast  Element 


2 


3.  Outer  Magnetosphere  and  Magnetopause  Forecast  Element 


•  determination  of  magnetic  field  in  outer  magnetospheric 
region 

•  prediction  of  location  of  approximate  pressure-bal anced 
magnetopause  surface 

•  determination  of  input  parameters  for  Hagnetosheath  and 
Bom  Shock  Forecast  Element 


4.  Hagnetosheath  and  Bom  Shock  Forecast  Element 

•  prediction  of  location  of  boM  shock 

•  prediction  of  plasma  and  magnetic  field  throughout 
magnetosheath 

•  determination  of  input  parameters  for  Forecast 

Presentation  Output  Element 

•  input  for  NSM  Model  and  other  components  of  forecast 

system 

5.  Forecast  Presentation  Output  Element 

•  numerical  values 

•  graphical  displays 


All  of  these  elements  Mere  successfully  embodied  in  a  preliminary 
SWT  model  that  Mas  developed  for  the  purpose  of  providing  utility  in  the 
interim  period  before  the  final  operational  SWT  model  became  available. 
In  particular,  the  preliminary  SWT  model  Mas  specifically  constructed  to 
provide  all  the  necessary  information  regarding  model  input  and  output 
that  Mould  be  required  for  real-time  operational  use  Mith  the  final 
operational  SWT  model.  The  preliminary  SWT  model  is  described  in  detail 
i n  Ref .  1 . 

Because  of  the  shortness  of  time  from  contract  start  to  delivery  of 
the  interim  SWT  model,  an  essential  element  that  Mas  not  possible  to  be 
carried  out  Mith  the  preliminary  SWT  model  Mas  the  detailed  evaluation 
of  SWT  model  forecasts  as  Mell  as  the  actual  SWT  model  validation 
employing  actual  satellite  data.  This  task  Mas  planned  to  be  the  focus 
of  the  next  technical  effort  to  be  carried  out  on  the  contract,  and  Mas 
conceived  to  consist  of  the  folloMing  subtasks: 
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•  establishment  of  a  attaints  data  base  that  real Istlcally 
simulates  real-time  operational  data  and  that  would  be 
suitably  for  SWT  nodal  forecast  testing 

•  SWT  Hodel  validation  by  detailed  comparison  of  nodal 
forecasts  with  this  observational  data  base 

•  Identification  of  elenents  In  the  solar  wind  transport 
nagnetosheath  forecast  systen  that  require  further 
snhancenent 

•  developnent  and  Incorporation  of  all  Inprovenents 
Identified  as  necessary  Into  the  final  operational  SWT 
nodel 


The  work  accomplished  on  these  subtasks  Is  described  below. 


2.  SATELLITE  DATA  BASE  FOR  VALIDATION  OF  SWT  FORECAST  NODEL 


A  study  was  carried  out  to  Identify  appropriate  candidate  archival 
satellite  data  sets  that  would  meet  the  requirements  for  use  In  the 
validation  program  for  the  SWT  model.  These  data  sets  would  necessarily 
be  composed  of  simultaneous  observations  from  at  least  two  spacecraft 
and  provide  concurrent  observations  of  both  oncoming  solar  wind  plasma 
and  magnetic  field  conditions  and  corresponding  conditions  In  the 
terrestrial  magnetosheath. 

As  a  result  of  this  Investigation.  It  was  found  that  the  most 
suitable  and  available  data  sets  that  meet  these  requirements  are  the 
1978*80  ISEE  2  and  3  observations.  A  large  library  of  archival  results 
from  these  ISEE  spacecraft  Is  currently  available.  However,  the  actual 
assembly  of  the  satellite  data  base  sufficient  for  validation  of  the  SWT 
model,  which  would  necessarily  Involve  the  correlation  of  a  large  number 
(of  the  order  of  *200)  of  simultaneous  ISEE  2  and  ISEE  3  data  sets  from 
this  library,  represents  a  substantial  effort.  This  Is  so  since  the 
Individual  ISEE  spacecraft  data  sets  have  to  be  searched  separately  to 
ascertain  data  continuity,  completeness  and  appropriateness,  and  then 
these  same  ISEE  data  sets  need  to  be  correlated  to  Insure  simultaneous 
coverage  of  both  solar  wind  and  magnetosheath.  Furthermore,  such  an 
effort  Is  peripheral  to  the  SWT  model  developments;  and  improvements 
needed  in  the  validation  effort.  Consequently,  in  order  to  ameliorate 
the  Impact  on  contract  resources,  we  have  established  a  joint 
collaborative  program  with  UCLA  and  Prof.  C.  T.  Russell  and  Dr.  6.  Le  of 
the  Geophysics  Dept .  Under  this  program,  RNA  and  UCLA  will  jointly 
identify  the  simultaneous  ISEE  data  sets.  UCLA  will  provide  the  actual 
numerical  ISEE  data  sets  to  RHA,  and  RHA  will  share  the  comparative 
results  with  UCLA.  The  result  of  this  collaborative  program  is  that  RHA 
will  acquire  a  sufficiently  extensive  satellite  data  base  to  validate 
the  SWT  model  at  a  minimum  expense  to  the  project. 

We  have  now  carried  out  the  Identification  of  184  simultaneous  ISEE 
solar  wind  and  magnetosheath  data  sets  that  are  both  appropriate  and 
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sufficient  for  the  SWT  model  validation  program.  These  data  sets 
Involve  ISEE  2  observations  of  the  magnetosheath  with  simultaneous  ISEE 
3  observations  of  the  oncoming  solar  wind.  Each  of  the  ISEE  2  data  sets 
Involves  at  least  one  magnetopause  crossing  and  may  contain  In  addition 
one  or  more  bow  shock  crossings.  The  time  Intervals  spanned  by  the  ISEE 
2  magnetosheath  traversals  vary  from  about  4  to  42  hours  and  the 
traversals  are  located  In  the  dayslde  magnetosheath.  For  these 
observations.  ISEE  3  Is  nominally  located  In  the  solar  wind  near  the  Li 
llbratlon  point.  All  the  observational  data  for  both  ISEE  3  and  ISEE  2 
for  both  plasma  and  magnetic  field  are  available  at  1  min.  Intervals. 

Table  1  provides  an  overview  of  the  combined  ISEE  3  and  ISEE  2  data 
sets.  A  summary  Is  provided  of  conditions  both  In  the  solar  wind  and  at 
the  Earth  (Dst,  dipole  tilt)  at  the  time  of  the  Initial  ISEE  2 
magnetopause  crossing.  The  data  sets  are  ordered  by  Increasing 
magnetosonic  Nach  number  and  span  the  range  from  2.1  £  Mms  £  8.2. 

Figures  1  and  2  below  provide  illustrations  of  sample  simultaneous 
data  sets  from  this  data  base.  Figure  2  displays  data  set  #003  which 
corresponds  to  ISEE  2  Orbit  314  inbound.  Hagnetosheath  data  from  ISEE  2 
are  shown  by  the  thin  lines,  while  corresponding  solar  wind  data  from 
ISEE  3  Is  superimposed  and  shown  as  thick  lines.  The  ISEE  3  data  for 
this  illustration  has  been  uniformly  time  shifted  forward  by  50  minutes 
to  allow  for  the  arrival  time  lag.  The  data  illustrated  in  Figure  1  are 
for  moderately  unsteady  solar  wind  plasma  conditions  as  evidenced  by  the 
observed  density  variations,  but  for  fairly  steady  INF.  These  data  are 
for  relatively  low  solar  wind  magnetosonic  Hach  number  conditions,  i.e. 
Mbs  -2.  Figure  2  displays  corresponding  ISEE  3  and  ISEE  2  data  for 

more  typical  solar  wind  magnetosonic  Hach  number  conditions,  i.e.  Mms  *» 
6.  with  conditions  in  the  solar  wind  fairly  steady  for  both  the  plasma 
and  IMF. 

The  current  status  of  this  ISEE  2  and  ISEE  3  data  base  is  that  the 
data  base  is  currently  archived  in  total  at  UCLA.  Selected  sets  of  the 
data  base  which  are  described  below  have  been  transferred  to  RHA  for 
initial  comparative  testing  of  the  SWT  forecast  model.  Final  transfer 
of  the  complete  archival  data  base  to  RHA  has  not  been  carried  out. 


3.  SWT  HODEL  VALIDATION  PROGRAH 

3.1  Previous  Validation  of  Elements  of  the  SWT  model 

The  primary  purpose  of  the  SWT  model  validation  program  is  to 
assess  and  qualify  the  accuracy  of  the  SWT  model  forecasts  in  situations 
that  simulate  as  closely  as  possible  the  real-time  operational 
environment  that  the  model  will  be  used  in.  Previous  studies  have 
served  to  successfully  verify  a  number  of  elements  embodied  in  the 
current  SWT  model.  For  example,  validation  of  the  methodology  employed 
in  the  development  of  the  approximate  pressure-balanced  magnetopause 
shape,  as  shown  in  Figure  3,  which  underlies  the  entire  basis  of  the 
gasdynamic  convected  field  model  employed  in  the  SWT  model  to  determine 
the  magnetosheath  plasma  and  magnetic  field,  was  made  early  on  in  Ref. 
2.  Similarly,  validation  of  the  accuracy  of  the  Briggs-Spreiter 
magnetopause  model  incorporated  in  the  SWT  model  was  successfully  made 
in  Ref.  3  by  comparison  with  other  theoretical  magnetopause  model 
results.  These  results  are  summarized  in  Figure  4.  An  important  early 
evaluation  of  the  gasdynamic  plasma  flow  model  employed  in  the  SWT  model 
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was  reported  in  Ref.  4  where  a  comparison  of  model -predicted  with 
laboratory-observed  bow  shock  shape  about  a  model  magnetopause  was  made. 
These  results  are  shown  in  Figure  5.  One  of  the  initial  evaluations  of 
detailed  plasma  properties  in  the  magnetosheath  as  predicted  by  the 
gasdynamic  model  was  made  by  comparison  with  observations  of  Pioneer  VI 
plasma  data  (Ref.  4)  as  shown  in  Figure  6.  and  also  with  IHP-1  data 
(Ref.  4)  and  is  shown  in  Figure  7.  Another  initial  validation  of  the 
gasdynamic  flow  model  but  this  time  using  magnetometer  rather  than 
plasma  observations  of  the  bow  shock  by  five  IMP  spacecraft  from  1963-68 
was  made  in  Ref.  5  and  is  shown  in  Figure  6.  A  more  recent  validation 
of  the  convected  magnetic  field  model  employed  in  the  SWT  model  was 
reported  in  Ref.  6  and  is  illustrated  in  Figure  9  where  a  comparison  of 
the  predicted  magnetic  field  along  the  Pioneer- Venus  spacecraft 
trajectory  in  the  Venusian  magnetosheath  was  made  with  PVO  magnetometer 
observations.  Another  successful  recent  validation  of  the  convected 
magnetic  field  model  is  reported  in  Ref.  7  and  illustrated  in  Figure  10. 
where  comparison  of  the  predicted  magnetic  field  along  the  Voyager  2 
spacecraft  trajectory  was  made  with  Voyager  2  magnetometer  observations. 
A  recent  summary  of  comparisons  of  various  elements  embodied  in  the  SWT 
model  with  observations  from  a  variety  of  spacecraft  and  solar 
wind/ pi  a net ary  interactions  throughout  the  solar  system  are  provided  in 
Refs.  7  and  8.  All  of  these  previous  comparisons  have  served  to  provide 
a  firm  basis  for  the  preliminary  successful  validation  of  the  accuracy 
of  a  number  of  elements  of  the  SWT  model. 


3.2  Current  Validation  Results  of  the  SWT  Model 

The  one  important  aspect  that  was  not  possible  to  investigate  in 
these  previous  studies  due  to  the  absence  of  an  appropriate  satellite 
data  base  was  the  use  of  continuous  oncoming  solar  wind  plasma  and 
magnetic  field  Information  as  input  to  the  SWT  model.  Prior  to  the 
current  study,  this  was  the  one  crucial  feature  of  the  operational 
forecasting  problem  that  had  never  been  previously  tested  against  the 
various  elements  embodied  in  the  present  SWT  model.  The  key  advantage  of 
the  SWT  model  as  employed  in  an  operational  mode,  however,  lies  in  this 
very  aspect.  That  is.  the  ability  of  the  currently- configured  SWT  model 
to  employ  short  time-interval  (of  the  order  of  1  min.)  solar  wind  input 
data  to  rapidly  predict  (in  a  few  CPU  secs.)  the  global  magnetosheath 
plasma  and  field  properties,  including  in  particular  the  location  of  the 
magnetopause  and  bow  shock  surfaces  and  conditions  at  specific 
spacecraft  orbital  locations  within  the  global  magnetosheath,  makes  the 
model  uniquely  valuable  for  real-time  operational  use. 

The  general  plan  for  the  SWT  model  validation  program  is  to  employ 
the  simultaneous  ISEE  3  solar  wind  data  as  input  to  the  SWT  model  and 
then  to  make  comparisons  of  the  plasma  and  magnetic  field  properties 
forecast  by  the  model  at  the  ISEE  2  spacecraft  orbital  locations  in  the 
magnetosheath  with  the  ISEE  2  observations  of  plasma  and  magnetic  field. 
In  order  to  qualify  the  SWT  model  for  operational  use,  it  will  be 
necessary  to  employ  approximately  half  of  the  184  simultaneous  ISEE  3 
and  ISEE  2  data  sets  contained  in  the  oata  base.  It  was  planned  that 
the  Initial  comparative  testing  would  be  made  for  selected  data  sets 
corresponding  to  solar  wind  conditions  believed  favorable  for  good  SWT 
model  forecasting.  Subsequent  comparisons  would  then  be  made  for 
selected  data  sets  corresponding  to  solar  wind  conditions  believed  to  be 
Increasing  less  favorable  for  good  forecasting,  i.e.  low  Mach  number, 
large  IMF,  very  unsteady  plasma  and  magnetic  field,  etc. 
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To  date,  we  have  now  carried  out  detailed  comparisons  of  the  SWT 
model  forecasts  with  the  ISEE  observations  for  6  of  these  simultaneous 
data  sets.  These  validation  studies  demonstrate  that  the  Phase  I  SWT 
model  forecasting  results  are  generally  in  good  agreement  with  the  ISEE 
2  magnetosheath  observations.  They  have  also,  however,  disclosed  a 
number  of  deficiencies. 

In  particular,  the  global  magnetic  field  forecast  of  the  Phase  I 
SWT  model  in  the  vicinity  of  the  magnetopause,  particularly  at  low  solar 
wind  Alfven  Mach  numbers,  was  found  to  display  a  consistent  systematic 
discrepancy  with  the  ISEE  observations.  The  source  of  a  major  part  of 
this  discrepancy  has  been  identified  as  due  to  a  previously- known 
spurious  magnetic  singularity  at  the  magnetopause  that  is  contained  in 
the  basic  gasdynamic  convected  magnetic  field  model.  A  provisional 
procedure  has  been  preliminarily  developed  to  diminish  the  effects  of 
this  singularity.  This  procedure  is  based  on  the  known  physical 
constraint  that  the  magnetosheath  magnetic  field  at  the  magnetopause 
cannot  exceed  the  local  confined  magnetospheric  magnetic  field  on  the 
interior  surface  of  the  magnetopause  at  that  location.  This  procedure 
has  been  initially  tested  against  the  ISEE  2  data,  and  the  forecasts  of 
the  magnetic  field  provided  by  this  modification,  as  shown  in  Figures  12 
to  23  below,  provide  a  significant  improvement  over  those  provided  by 
the  Phase  I  SWT  model . 

In  addition.  Initial  comparisons  indicate  that  superior  SWT  model 
forecasts  are  achieved  at  low  solar  wind  Alfven  Mach  numbers  by 
employing  the  magnetosonic  rather  than  sonic  Mach  number  as  input  to  the 
SWT  model.  Additionally,  it  has  been  found  that  employing  the 
magnetosonic  Mach  number  does  not  alter  the  agreement  of  the  SWT  model 
forecasts  for  situations  at  high  solar  wind  Alfven  Mach  numbers.  In  all 
of  the  comparisons  shown  below  in  Figures  11-23,  the  magnetosonic  Mach 
number  as  been  employed  as  input  to  the  gasdynamic  calculation  in  the 
SWT  model . 

In  Figures  11  and  12.  we  provide  results  of  the  first  comparisons 
of  the  SWT  model  forecasts  and  ISEE  2  observations  for  ISEE  2  orbit  145 
outbound.  This  data  set  (#178)  was  selected  as  one  of  the  initial 
candidates  for  the  validation  study  as  it  corresponded  to  solar  wind 
conditions  that  were  believed  favorable  for  good  SWT  model  forecasts; 
that  is,  high  magnetosonic  Mach  number  (Mms  >7.1),  B*  northward,  steady 

solar  wind  plasma  conditions,  fairly  steady  IMF,  and  single  bow  shock 
and  magnetopause  crossings  by  ISEE  2.  The  plot  on  the  left  of  Figure  11 
provides  a  view  of  the  ISEE  2  spacecraft  orbit  as  seen  in  the  solar  wind 
(X.R)  coordinate  system  fixed  to  the  magnetopause  that  is  employed  in 
the  global  magnetosheath  computations.  Such  a  view  of  the  spacecraft 
trajectory  Is  quite  informative  as  it  immediately  exhibits  the  magnitude 
of  the  effects  caused  by  unsteadiness  in  the  solar  wind  dynamic  pressure 
and  solar  wind  direction.  These  are  revealed,  respectively,  in  this 
view  of  the  spacecraft  trajectory  by:  (1)  abrupt  radial  jumps  in  the 
spacecraft  trajectory  corresponding  to  the  sudden  inflation  or  deflation 
of  the  magnetopause  obstacle  size  in  response  to  solar  wind  dynamic 
pressure  changes,  and  (2)  abrupt  lateral  excursions  of  the  spacecraft 
trajectory  corresponding  to  sudden  changes  in  the  solar  wind  direction. 
As  is  evident  in  the  plot  in  Figure  11.  for  the  time  that  the  ISEE  2 
spacecraft  spent  in  the  magnetosheath  during  this  orbit,  the  direction 
of  the  solar  wind  underwent  several  abrupt  changes.  The  comparison  of 
the  prediction  from  the  SWT  model  for  the  velocity  component  Vx  in  GSE 


7 


coordinates  with  the  ISEE  2  fast  plasma  analyzer  displays  very  good 
agreement  for  the  location  of  both  the  bow  shock  and  magnetopause 
surfaces,  and  very  good  agreement  for  the  magnitude  of  the  velocity 
component  for  most  of  the  magnetosheath  crossing.  Some  disagreement 
does  appear  between  SWT  predictions  and  observations  in  the  vicinity  of 
the  magnetopause  surface.  We  have  no  explanation  for  this.  As  will  be 
seen  from  the  comparisons  for  the  other  five  case  studies,  the 
discrepancies  of  SWT  model  forecasts  with  the  ISEE  2  fast  plasma 
analyzer  observations  for  velocity  exhibit  no  systematic  trends.  In 
some  cases,  the  agreement  is  essentially  perfect,  while  for  others  some 
discrepancies  may  appear.  These  discrepancies  might  in  fact  be 
attributable  primarily  to  the  instrument  since  we  have  discovered  what 
appears  to  be  a  systematic  error  in  the  magnitude  (underreported  by 
almost  a  factor  2)  of  the  plasma  density  reported  by  the  ISEE  2  fast 
plasma  analyzer  for  all  of  the  case  studies  that  have  been  undertaken. 
Figure  12  provides  corresponding  comparison  of  the  SWT  model  forecast 
and  the  ISEE  2  magnetometer  for  the  magnetic  field  components  in  GSE 
coordinates.  As  with  the  plasma  velocity  comparison,  the  magnetic  field 
comparisons  are  very  good  for  most  of  the  magnetosheath  region,  with 
some  disagreement  in  the  vicinity  of  the  magnetopause.  These 
comparisons  confirm  the  good  agreement  of  the  locations  of  both  the  bow 
shock  and  magnetopause  surfaces  already  observed  in  the  plasma  velocity 
comparisons  in  Figure  11.  To  provide  some  Indication  of  the  steadiness 
of  the  solar  wind  parameters  during  the  ISEE  2  magnetosheath  crossing 
for  this  orbit,  we  have  prepared  Figure  13  which  provides  the  variation 
of  the  solar  wind  sonic,  Alfven  and  magnetosonic  Hach  numbers  during  the 
ISEE  2  magnetosheath  traversal.  Also  shown  in  Figure  13  is  the 
pressure-balanced  magnetopause  nose  radius  Rp.  which  provides  an 
indication  of  the  inflation  and  deflation  undergone  by  the  magnetopause 
during  this  time  Interval.  As  can  be  observed  from  Figure  13.  while  the 
solar  wind  conditions  are  reasonably  steady  during  this  time  interval, 
they  are  by  no  means  absolutely  quiescent. 

Figures  14  and  15  provide  corresponding  comparisons  of  SWT  model 
forecasts  and  ISEE  2  observations  for  orbit  178  inbound.  This  data  set 
(#181)  was  selected  as  one  of  the  candidates  for  the  initial  validation 
study  since  although  this  orbit  contains  solar  wind  conditions  that  were 
considered  mostly  favorable  for  good  SWT  model  forecasts,  it  also 
contains  some  solar  wind  conditions  believed  unfavorable  that  would 
subsequently  serve  as  a  further  test  of  the  SWT  model.  For  this  case, 
the  solar  wind  conditions  are:  high  magnetosonic  Hach  number  (M^  « 

8.0),  Bz  northward,  moderately  unsteady  solar  wind  plasma  conditions, 
moderately  unsteady  IHF,  and  multiple  bow  shock  crossings  were  observed 
by  ISEE  2.  In  Figure  14,  we  observe  from  the  plot  of  the  ISEE  2  orbit 
In  solar  wind  coordinates  that  due  to  inflation  of  the  magnetopause  the 
spacecraft  essentially  loiters  about  the  bow  shock  resulting  in  the 
multiple  bow  shock  crossings  observed  in  the  plasma  velocity  comparisons 
shown  In  the  plot  on  the  right.  The  ISEE  3  data  gap  during  the  period 
when  the  spacecraft  initially  crosses  the  bow  shock  undoubtedly  causes 
the  discrepancy  for  the  first  predicted  bow  shock  crossing.  However, 
the  prediction  of  the  later  multiple  bow  crossings  by  the  SWT  model  and 
the  magnitude  of  the  velocity  component  throughout  the  magnetosheath  is 
In  excellent  agreement  with  the  ISEE  2  fast  plasma  analyzer 
observations.  Figure  15  also  confirms  this  high  level  of  agreement  for 
the  magnetic  field  components  as  well,  and  In  addition  shows  that  even  a 
difficult  detail,  such  as  the  multiple  bow  shock  crossing  at  about  14.0 
hours  UT.  is  remarkably  captured  by  the  SWT  model  magnetic  field 
forecasts.  Both  the  level  and  variation  of  the  all  of  the  magnetic 
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field  components  are  exceptionally  well  predicted  for  this  case  In  which 
there  exists  moderate  unsteadiness  in  both  the  oncoming  solar  wind 
plasma  flow  as  well  as  the  IHF. 

Figures  16  and  17  display  the  analogous  comparisons  for  ISEE  2 
orbit  139  outbound.  This  data  set  (#145)  was  selected  as  one  of  the 
candidates  for  this  initial  validation  study  as  it  corresponded  to  solar 
wind  conditions  that  were  also  believed  favorable  for  good  SWT  model 
forecasts:  that  is.  high  magnetosonic  Hach  number  (Mms  -  5.6).  B2 

northward,  steady  solar  wind  plasma  conditions,  fairly  steady  IHF,  and 
single  bow  shock  and  magnetopause  crossings  by  ISEE  2.  In  Figure  16,  we 
note  from  the  plot  of  the  ISEE  2  trajectory  in  solar  wind  coordinates 
that  a  high  number  of  directional  changes  in  the  oncoming  solar  wind 
occur  during  the  ISEE  2  magnetosheath  traverse  of  this  orbit. 
Nevertheless,  the  comparison  of  the  SWT  model  forecast  and  ISEE  2 
observed  axial  velocity  component  exhibits  very  good  agreement  for 
location  of  both  bow  shock  magnetopause  crossings.  In  addition,  the 
trend  of  the  variation  of  this  velocity  component  is  well  predicted, 
although  the  magnitude  appears  to  be  off  by  a  constant  factor  throughout 
the  magnetosheath.  As  before,  we  have  no  explanation  for  this 
discrepancy.  In  Figure  17,  however,  the  three  components  of  predicted 
magnetic  field  display  excellent  agreement  with  the  ISEE  2  observations 
across  the  entire  magnetosheath.  This  is  a  key  reason  why  we  suspect 
errors  in  the  ISEE  2  fast  plasma  analyzer  data.  Recall  that  in  the  SWT 
model,  the  convected  magnetic  field  is  the  last  global  quantity  to  be 
calculated.  This  is  so  since  the  convected  field  model  incorporated  in 
the  SWT  forecast  model  requires  knowledge  of  both  the  velocity  and 
density  fields  throughout  the  magnetosheath  before  the  computation  of 
the  magnetic  field  can  be  initiated.  Consequently,  any  error  in  the 
flow  field  computation  or  any  inadequacy  in  either  the  flow  field  or 
magnetic  field  computation  would  be  displayed  most  prominently  in  the 
magnetic  field  comparisons.  However,  for  this  case,  as  well  as  others 
studied  here  in  which  a  discrepancy  appears  in  the  comparisons  between 
the  SWT  model  predictions  and  ISEE  2  observations  for  the  magnetosheath 
velocity,  this  discrepancy  in  velocity  is  not  carried  over  into  the 
magnetic  field  comparisons  which  almost  uniformly  display  excellent 
agreement.  We  note,  furthermore,  that  the  agreement  between  SWT  model 
forecast  and  observation  is  also  quite  good  in  all  cases  studied  for  the 
locations  of  both  the  bow  shock  and  magnetopause  surfaces.  It  is  not 
possible  for  the  SWT  model  forecasts  to  be  accurate  in  predicting  the 
locations  of  the  bounding  surfaces  of  the  magnetosheath  as  well  as  to  be 
accurate  in  the  predictions  of  the  vector  magnetic  field  throughout  the 
magnetosheath  while  simultaneously  containing  a  significant  error  in  the 
magnetosheath  vector  velocity.  Consequently,  we  believe  that  there 
exist  errors  in  the  ISEE  2  plasma  analyzer  data  for  the  vector  velocity 
as  well  as  for  the  density.  Further  study  of  this  aspect  of  the  ISEE  2 
observations  is  clearly  needed. 

In  Figures  18  thru  23.  analogous  comparisons  are  made  between  the 
SWT  model  forecasts  and  ISEE  2  magnetosheath  observations  for  three 
additional  case  studies  all  of  which  involve  low  solar  wind  magnetosonic 
Mach  number,  i.e.  2.2  £  M^  s  3.0.  These  data  sets  (#009, #007 .#003) 
were  purposely  selected  as  candidates  for  this  initial  validation  study 
in  order  to  provide  a  further  challenge  to  the  SWT  model  since  the 
assumptions  under  which  the  basic  gasdynamic  convected  magnetic  field 
computational  model  used  to  determine  the  global  magnetosheath  plasma 
and  magnetic  field  properties  in  the  SWT  model  is  based  become  highly 
strained  at  low  magnetosonic  Mach  number.  However,  with  the  use  of  the 
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solar  wind  magnetosonlc  rather  than  sonic  Hach  number  In  the  global 
gasdynamlc  plasma  calculation,  the  predicted  results  for  both  plasma  and 
magnetic  field  properties  In  the  magnetosheath  for  all  3  of  these  case 
studies  are  In  very  good  agreement  with  the  ISEE  2  observations.  This 
result  was  not  anticipated,  and  warrants  further  study  to  evaluate  the 
reasons  underlying  this  agreement.  Figures  18  and  19,  which  correspond 
to  ISEE  2  orbit  319  outbound,  display  results  for  the  longest  T'PE  2 
magnetosheath  traversal  Interval  studied  In  all  of  the  case  es 
undertaken,  i.e.  41  hours.  Figure  18  shows  very  good  agreement  sen 
the  SWT  model  forecast  and  ISEE  2  observation  for  the  velocity 
component,  and  Figure  19  confirms  that  agreement  for  the  vector  magnetic 
field.  Note  In  particular  the  extreme  variation  In  the  Bz  component 
displayed  In  the  ISEE  2  data  In  Figure  19  and  which  Is  accurately 
captured  by  the  SWT  model  forecasts.  Figures  20  and  21  display 
analogous  comparisons  for  ISEE  2  orbit  331  Inbound.  Figure  20  shows  a 
remarkable  agreement  between  the  SWT  model  forecasts  and  ISEE  ? 
observations  for  the  axial  velocity  component,  with  the  only  exception 
being  the  SWT  model  forecast  of  a  bow  shock  crossing  at  about  7:30  UT 
which  does  not  appear  In  the  observations.  Figure  21  which  displays  the 
comparison  of  the  vector  magnetic  field  components  exhibits  the  same 
excellent  agreement  between  model  forecasts  and  observations.  Figures 
22  and  23  display  the  final  comparisons,  and  correspond  to  solar  wind 
conditions  for  the  lowest  magnetosonlc  Kach  numbers  of  all  the  case 
studies  undertaken.  Again,  the  agreement  between  the  SWT  model 
forecasts  and  the  ISEE  2  observations  for  the  axial  velocity  component 
shown  In  Figure  22  Is  very  good,  as  Is  that  for  the  vector  magnetic 
field  Illustrated  In  Figure  23.  The  only  feature  of  these  three  case 
studies  that  displays  any  systematic  discrepancy  between  model  forecast 
and  observation  at  all  Is  that  the  magnitude  of  the  predicted 
magnetosheath  magnetic  field  at  low  magnetosonlc  Hach  numbers  Is 
somewhat  underpredicted  compared  to  the  observations.  A  correction 
parameterized  by  magnetosonlc  Hach  number  and  based  on  an  extended 
series  of  cases  studies  should  be  developed  In  order  to  Improve  this. 


4.  CONCLUSIONS  AND  RECOHHENDATIONS 
4.1  Status  of  SWT  Validation  Program 

The  observational  data  base  required  for  the  validation  of  the  SWT 
forecasting  model  has  now  been  established.  The  data  base  Is  comprised 
of  a  total  of  184  ISEE  2  and  simultaneous  ISEE  3  data  sets  that  Involve 
solar  wind  plasma  and  field  measurements  from  the  ISEE  3  spacecraft  and 
corresponding  plasma  and  field  magnetosheath  measurements  from  the  ISEE 
2  spacecraft.  Case  studies  to  Investigate  the  accu  acy  of  the  SWT  model 
forecasts  have  been  Initiated  by  using  selected  sets  from  this  data 
base.  Comparative  results  have  currently  been  obtained  for  6  of  these 
ISEE  2  and  ISEE  3  data  sets.  These  selected  sets  have  Involved  solar 
wind  conditions  anticipated  to  be  both  favorable  and  unfavorable  for 
good  SWT  forecasts.  The  preliminary  results  from  these  case  studies 
have  In  general  been  very  favorable.  Even  from  this  small  number  of 
cases  studied  so  far,  however,  several  Important  features  in  both  the 
plasma  and  magnetic  field  procedures  employed  In  the  Interim  Phase  I  SWT 
model  have  been  Identified  that  require  Improvement.  Provisional 
procedures  have  been  developed  as  far  as  possible  during  this  period  to 
Improve  these  methods,  but  require  both  further  development  and  study. 
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4.2  Final  Qualification  of  SHT  Forecasting  Model 

At  the  present  development  stage,  the  SHT  forecasting  model  remains 
unqualified  for  operational  use.  The  study  to  establish  figures  of  merit 
for  the  model  that  would  provide  estimates  of  the  quality  of  the  real¬ 
time  SHT  forecasts  in  an  operational  environment  has  not  been  initiated. 

The  primary  technical  tasks  remaining  to  complete  the  SHT  model 
validation  program  and  qualify  the  model  for  operational  use  consists  of 
the  following: 

-  completion  of  the  SHT  model  comparative  forecasts  for  a  minimum 
of  half  of  the  currently  established  ISEE  2  and  ISEE  3 
observational  data  base 

*  development  and  incorporation  of  all  improvements  identified  as 
necessary  from  the  above  comparative  study  for  the  plasma  and 
magnetic  field  computational  models  embodied  in  the  Phase  I  SHT 
model 

-  establishment  of  figures  of  merit  for  the  SHT  model  forecasts  as 
a  function  of  solar  wind  parameters 
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Table  1.  Overview  of  Combined  ISEE  Z,  ISEE  3  Data  Sats  For  Phase  I 
SWT  Nodal  Validation  Program:184  ISEE  Z  Magnetosheath 
Traverses  With  MP/8ow  Shock  Crossings ;0rdered  by 
Magnetosonlc  Mach  Nuaiber 


nrmcncnrerigTiEicaag 


Tabic 


Data  Sc*  #003:  ISEE  3  Solar  Wind  +  ISEE  2  Magnetosheath  Data 


Figure  1.  Illustration  of  Sample  Simultaneous  ISEE2  and  ISEE  3  Data 
Sat  For  ISEE  2  Orbit  314  Inbound 
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Data  Set  #145:  ISEE  3  Solar  Wind  +  ISEE  2  Magnetosheath  Data 

ISEE  2  (Thin  line)  and  ISEE  3  (Thick  line).  Delta  Time  =  50  min 


BX  GSE 
nT 


BY  GSE 
nT 


BZGSE 
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DEN 
cm- 3 
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Figure  2.  Illustration  of  Sample  Simultaneous  ISEE2  and  ISEE  3  Data 
Set  For  ISEE  2  Orbit  139  Outbound 


I  0  -I 


X/D 

Comparison  of  exact  and  approximate  pressure  distributions  on  magnetosphere  boundary 


Figure  3.  Early  Validation  of  Elements  of  SWT  Model:  Evaluation  of 
Approximate  Newtonian  Pressure  Formula  by  Comparison  With 
Gasdynamlc  Flow  Computations 
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Magnetopause  Model  by  Comparison  With  Other  Theoretical 
Magnetopause  Results 


Figure 
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Comparison  of  calculated  and  obaerued  portions  of  bow  wave  of  modal  magnetosphere 
in  flighl  at  Mach  number  4.69  through  argon  (7  —  3/3).  Equatorial  plana 
(Sruma  and  Hyttt,  1963),  pjp*  —  .1. 


Early  Validation  of  Elements  of  SHT  Model:  Evaluation  of 
Gasdynamlc  Flow  Model  by  Comparison  With  Laboratory 
Experiment 


19 


20 


Figure 


THEORY 


IMP-I  plasma  probe  measurements  of  shock  wave  and  magnetosphere  boundary  crossings 


7.  Early  Validation  of  Elements  of  SWT  Model:  Evaluation  of 
Gasdynamlc  Flow  Model  by  Comparison  With  Observations  of 
IMP-1 
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VMM. 

Point* 


Position  of  th«  bow  ihock  in  the  wlar  ecliptic  plan*  as  determined  bp  meaewo* 
ta  ob  live  Imp  apacecraft,  963-1968.  Crown  represent  the  areraf*  loc«tioa  oo  iadividaal 
tad  the  tolid  line  hyperbole  represent*  the  beet-fit  curve  to  the  point*.  jZaa|  <  7  Jl» 
t*  hev*  been  rotated  by  4*  to  remove  the  effects  of  aberration  duo  to  the  earth 'i  motion 


about  the  ana.  The  tin* 


to  beyond  the  av 


*f  neneneWy  distant  bow  ancafc  leeolione. 


Figure  0.  Early  Validation  of  Claaenta  of  SHT  Nodal:  Evaluation  of 
Gasdynamlc  Flow  Modal  by  Comparison  With  Hagnetometar 
Observations  of  Bow  Shock  by  Five  IMP  Spacecraft ; 1963- 1968 


Figure  10.  More  Recent  Validation  of  Elements  of  SWT  Model:  Evaluation 
of  Convected  Magnetic  Field  Model  at  Extreme  Low  Density 
Flow  Situation  (n  -  0.005/cm3)  by  Comparison  With  Voyager  2 
Observations  at  Neptune 


Dad  Sat  #171 


ISEE  2  Orbit  145  Outbound  as  Viewed  froae  Solar  Wind 

(XJt)  Coordinate  Syatem  Fixed  to  bhfMlepMM  • 


Flgura  II.  Comparison  of  SWT  Modal  Forecasts  and  ISEE  2  Observations 
for  Orbit  145  Outbound : ISEE  2  Orbit  In  Solar  Wind 
Coordinates  and  Comparison  of  GSE  Velocity  Component  Vx 


Flgur#  12.  Coaparlson  of  SWT  Model  Forecasts  and  ISEE  2  Observations 
for  Orbit  145  Outbound:  GSE  Magnetic  Field  C opponents 
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Figure  13  Illustration  of  Variation  of  Solar  Wind  Sonic,  Alfvan  and 
Hagnatosonlc  Hach  Number  and  Magnetopause  Nose  Radius  During 
ISEE  2  Orbit  14S  Magnetosheath  Traversal 
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DataSet  till 

I  SEE  2  Orbit  171  Inbound  as  Viewed  from  Solar  Wind 
(XJR)  Coordinate  System  Fixed  to  Magnetopause 


Hflur*  14.  Coaparlson  of  SWT  Model  Forecasts  and  ISEE  2  Observations 
for  Orbit  178  Inbound:lSEE  2  Orbit  In  Solar  Hind  Coordinates 
and  Coaparlson  of  GSE  Velocity  Component  Vx 
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Flgur*  15.  Comparison  of  SWT  Model  Forecasts  and  ISEE  2  Observation* 
for  Orbit  170  Inbound:  GSE  Magnetic  Field  Components 
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Figure  IS.  Comparison  of  SWT  Model  Forecasts  and  ISEE  2  Observations 
for  Orbit  139  Outbound : ISEE  2  Orbit  In  Solar  Hind 
Coordinates  and  Comparison  of  GSE  Velocity  Component  Vx 
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Flgurt  17.  Comparison  of  SWT  Nodal  Forecasts  and  ISEE  2  Observations 
for  Orbit  139  Outbound:  GSE  Magnetic  Field  Components 
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Figure  18.  Comparison  of  $HT  Model  Forecasts  and  ISEE  Z  Observations 
for  Orbit  319  Outbound : ISEE  Z  Orbit  In  Solar  Hind 
Coordinates  and  Comparison  of  GSE  Velocity  Component  Vx 
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Figure  19.  Comparison  of  SMT  Modal  Foracasts  and  ISEE  2  Obsarvatlons 
for  Orbit  319  Outbound:  GSE  Magnetic  Field  Components 
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Data  Set  MT 


Modnwly  MUtatfr  mIv  Md  phono  cantfiwM 


BEE  2  Orbit  331  Inbound  ns  Viewed  from  Solar  Wind 
(XJR)  Coordinate  System  Fixed  to  Magnetopause 


_  Comparison  of  SWT  Model  and  ISEE  2  Orbit  331  Inbound  Data: 

VV**Z2  GSE  Velocity  Component  V* 


Figure  ZO.  Comparison  of  SWT  Model  Forecasts  and  ISEE  Z  Observations 
for  Orbit  331  Inbound: ISEE  2  Orbit  in  Solar  Hind  Coordinates 
and  Comparison  of  GSE  Velocity  Component  Vx 
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•  Moderately  unsteady  tolar  wind  plasma  conduons 

•  Moderately  atmeady  IMF  condnons 

•  Multiple  BS  end  ejagle  MF  croesum 


Flgur*  21.  Comparison  of  SHT  Model  Forecasts  and  ISEE  2  Observations 
for  Orbit  331  Inbound:  GSE  Magnetic  Field  Components 
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Data  Sot  0003 

1SEE  2  Orbit  314  Inbound  as  Viewed  from  Solar  Wind 
(X,R)  Coordinate  System  Fixed  to  Magnetopause 


ftwpdtfif  Klecttd  to  invotvt  tont  nkr  wind  co 
believed  unfavorable  for  good  SWT  model  foracam 

Ma(neto«on«c  Mach  number  low.  •  2  2 


Vyi*2i 


Comparison  of  SWT  Model  and  ISEE  2  Orbit  314  Inbound  Data: 
GSE  Velocity  Component  V* 


Range  of  Predicted  Bow 
Shock  Locations  From 
18:20  UT  to  19:50  UT/ 


36 


Flgur#  23.  Comparison  of  SWT  Model  Forecasts  and  ISEE  2  Observations 
for  Orbit  314  Inbound:  GSE  Magnetic  Field  Components 
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